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We propose a scheme to couple NV centers in diamond through coplanar waveguide resonators. 
The central conductor of the resonator is split into several pieces which are coupled strongly with 
each other via simple capacitive junctions or superconducting Josephson junctions. The NV centers 
are then put at the junctions. The discontinuity at the junctions induces a large local magnetic field, 
with which the NV centers are strongly coupled to the circuit resonator. The coupling strength g 
between the resonator and the NV center is of order of g/2n ~ 1-30 MHz. 



PACS numbers: 42.50.Pq, 03.65.-w, 03.67.-a, 37.30. +i 

A nitrogen-vacancy (NV) center in diamond consists 
of a nitrogen atom substituting a carbon atom and a va- 
cancy trapped adjacent to the substitutional nitrogen. In 
its ground state, the negatively charged NV center has 
a spin triplet, which is separated by optical transitions 
from the excited states. The long spin coherence time 
and fairly easy optical initialization and read out of the 
ground spin state make the NV center an excellent can- 
didate for quantum information processor and quantum 
information storage [1]. The coherent manipulation of 
the single spin or multiple spins of electrons and nuclei 
within a single NV center [2-4] or locally interacting NV 
centers [5] have been demonstrated experimentally. To 
build a scalable quantum information processors (QIP), 
however, controlled coupling between distant NV centers 
is yet to be achieved. 

A superconducting coplanar waveguide resonator, de- 
fined by a single centimeters-long central conductor be- 
tween two ground half-planes, has been successfully used 
as a quantum bus for charge-based superconducting 
qubits [6, 7] by exploiting a strong electric dipole cou- 
pling [8, 9]. Building a hybrid quantum device using the 
superconducting resonator as a quantum bus for the spin 
qubits (including NV centers) is desirable for the scalable 
QIP, because then we can take advantage of the scalabil- 
ity and low dissipation of the circuit-QED system as well 
as the long coherence time of spin qubits. However, it 
has been limited only to collective excitation of the spin 
ensemble due to a small magnetic dipole coupling of a 
single spin qubit with the resonator of an order of 10 
Hz [10-14]. Using a flux qubit as a mediator between the 
NV center and the resonator is suggested to enhance the 
coupling strength recently, but then the short coherence 
time of the flux qubit limits the coherence time of the 
entire system [15]. 

In this paper, we describe how to realize a quantum 
bus for distant NV centers using a series of superconduct- 



ing resonators coupled by either Josephson or capacitive 
junctions. We note that inserting a Josephson junction 
into the central conductor of the resonator enhances local 
magnetic field by a factor of 10 4 [16, 17]. Interestingly, 
a simple capacitive junction can also enhance the local 
magnetic field by the same factor as demonstrated be- 
low, which would allow us to circumvent the difficulties 
of fabricating many Josephson junctions in the central 
conductor. By putting single NV center at each junc- 
tion, the magnetic dipole coupling strength between the 
resonator and the NV centers reaches g/2ir ~ 1-30 MHz. 
The interaction between NV centers can be achieved by 
exchanging virtual photons and can be turned on and off 
by bringing them in and out of resonance with the res- 
onator frequency by the external magnetic field. Thus 
quantum gate operations can be performed in the same 
manner as the circuit-QED system [6, 7], An important 
difference is that the NV center qubit can be measured 
optically instead of being measured dispersively through 
the cavity. It means the high-loss cavity used for the fast 
measurement in the circuit-QED experiments [6, 7] is no 
longer necessary. We also note that our scheme offers 
orders of magnitude stronger coupling strength than re- 
cently proposed nanomechanical resonators-based quan- 
tum bus [18, 19] for the NV centers as well as schemes 
that use the flux qubits [15, 20], which enables us to 
perform fast quantum gate operations. Moreover, re- 
cent experimental demonstration of the coupling between 
an ensemble of NV centers and the superconducting res- 
onator [13] along with experiments realizing the super- 
conducting resonator with the Josephson junction in- 
serted in the central conductor [17, 21] indicate a fea- 
sibility of our scheme. 

Resonator. — We start with a general description of 
the resonator. We insert N junctions into the central 
conductor of the superconducting coplanar waveguide, 
essentially breaking the resonator into N + 1 subres- 
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FIG. 1: (color online) A schematic diagram of the coplanar 
waveguide resonator with (a) capactive or (b) Josephson junc- 
tions, (c) The energy-level structure of the NV center in dia- 
mon. 



onators. The adjacent subresonators are coupled via 
the junctions. The central conductor of the rth sub- 
resonator occupies the space X2 r < x < a^r+i with 
length £ r = X2r+i — %2r {r = 0, 1, 2, • • • , N) and has ca- 
pacitance Co per unit length. The junctions separating 
the adjacent subresonators have gap D r = X2 r — X2 r -i 
(r = 1,...,N) and junction capacitance Cj. We will 
assume that t r ~ 10 mm and D r ss 1-10 nm. 

The Lagrangian Cr governing the dynamics of the res- 
onator has two parts: Cr = Cq + Cj. The first part Cq 
describes the decoupled subresonators and is given by [8] 



Co 
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JV 



dx[(d t <f>) 2 -v 2 (d x <f>) 2 ] (1) 



where v ~ 10 8 m/s is the propagating velocity of the 
electromagnetic wave in the coplanar waveguide. Phys- 
ically, the field <p(x,t) is proportional to the magnetic 
flux and related to the local electric potential V(x,t) by 
dt4>(x,t) = V(x, t) The second part Cj comes from the 
coupling between adjacent subresonators: 



1 



TV 



(2) 



Here the variable <p r (i) is the magnetic flux within the 
junction and related to the electric potential difference 
V r (t) across the junction by d t <p r (t) = V r (t) . The first 
term in Eq. (2) is thus responsible for the electric energy 
stored in the junction and the second, the magnetic en- 
ergy. Up = y / 2EcEj/h w 27T x 10 GHz is the Josephson 
plasma frequency, where Ec = (2e) 2 /2Cj and Ej are the 
charging and Josephson coupling energy of the junction, 
respectively. The model (2) for Josephson junctions is 
valid only in the range ksT <C hu} p <C Ej whereas one 
can put uip = in (2) for capacitive junctions. Note 
that for large D r the junction capacity Cj (as well as the 
Josephson coupling Ej) becomes very small compared 
with £ r Co ~ lpF, and the coupling between adjacent 



subresonators are negligible except for a small red-shift 
of order oiCj/i T C [8]. For D r w 1-10 nm, Cj/£ r C ~ 1 
and the coupling is strong. 

The field <p(x, t) defined inside subresonators and tp r (t) 
defined across the junctions are not independent. They 
are related to each other by the current conservation 



v 2 C 



(d 2 +UJp)ip r (t) = d x (f)(x2r,t) = d x <f>(X2 r -l,t) (3) 



Following Ref. [16], we expand the field <p and <p r in the 
normal modes as 

oo oo 
<t>{x,t) = ^ 4>7n(t)^ m (x) , <p r (t) = ^ <Pm(t)^r,m- 



ro=Q 



m—0 



(4) 

The eigenfunctions tp m {x) should satisfy the time- 
independent Schrodinger equation 

(d 2 + k 2 m )^ m (x) = (x 2r <x< x 2r +i) (5) 

for some wave numbers k m . Let w m = vk m . The current 
conservation relation (3) now reads 

C 

-— -(Wp - u4jA rm = d x 1p m (x2r) = d x t/j m (x2r~l) (6) 

Without loss of generality, we choose the normalization 

N rX2 r+ i N 

^ / dx Tpm^n+Cj ^ A r , m A r ^ n = CsS mn , (7) 

r=0 Jx 2r r=1 

where Cs = J2r=o UCq + NCj is the total capacitance of 
the resonator. Putting the normal mode expansion (4) 
into Eqs. (1) and (2) and imposing the conditions (5), 
(6), and (7), one can rewrite Cr into the simple form 



(8) 



By introducing a momentum 9 m — Csdt4> m conjugate to 
4> m , we write the Hamiltonian of the resonator as 

(9) 

s / 

We then quantize it by the canonical commutation re- 
lation [4> m ,0 m ] = ih. It is customary to introduce the 
annihilation and creation operators of the normal modes 
by the relations 



<t>m = Y ^~^ ( a m + °™)' ( 10a ) 

m = i\ Z S ( a L~ a m) ; (10b) 



in terms of which the resonator Hamiltonian reads 



(11) 



m=0 



Due to the nonlinearity induced by the junctions, io n 
determined by (6) is not an integer multiple of the uq. 
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Coupling to NV centers. — We now describe the di- 
amond NV centers and their coupling to the resonator. 
The ground-state triplet (spin 1) of the NV center has 
a level splitting e/2ir = 2.88 GHz due to the spin-spin 
interaction [22, 23], and thus is described by the Hamil- 
tonian 



NV 



N 



(12) 



where SV,z is the spin z component of the NV center at 
the rth junction. 

The ground-state spin triplet is coupled magnetically 
to the resonator, which is governed by the coupling 
Hamiltonian 



N 

£ 

r=l 



9e^BS r ^ z B r 



(13) 



where B r = ip r /D r W (W is the distance between the 
central conductor and the ground plates) is the local 
magnetic field at the junction r, g e ~ —2 is the electron 
^-factor, and fis is the Bohr magneton. 

Using the annihilation and creation operator defined 
in Eq. (10), the coupling Hamiltonian is written as 



^ ^ 9m, r 

ra.r 



(a 



Here we have defined the coupling constant 

.9m, r 9ef^B B m _ r 



(14) 



(15) 



where B m r is the root-mean square value of the local 
magnetic field in the mode m at the junction r. By ad- 
justing one of the mode frequencies at resonance with 
the level splitting of the NV centers (cj m w e), one can 
selectively couple the mode m to the NV centers. 

Putting all together, the total Hamiltonian is given by 
the sum H = Hb + H^v + H g . The description above is 
completely general (in principle) for any number of junc- 
tions and NV centers. Equations (5) and (6) determines 
the possible spectrum uj m (m = 0, 1, 2, • • • ) of the coupled 
subresonators whereas the normalization (7) determines 
the magnitude of the magnetic field at the junctions. Be- 
low we demonstrate the cases with single NV center and 
two NV centers. 

Single NV center. — We assume that the junction is 
located at x = 0, and the two subresonators are identical 
with length £. The wave function il) m {x) of the normal 
model m has the form[10] 

y - cos[k m (x - t)\ [x > 0) . 

The wave numb er k m and t he frequency cj m — v k ni 
are determined by the current conservation relation (6), 
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FIG. 2: (color online) (a,b) Resonance frequencies ui m of the 
modes m = 0, 1, 2 (from bottom to top) as a function of the 
overall length £ of the resonator for uj p /2tt = 10 GHz (a) and 
ui p — (b). (c,d) The coupling strength g m to the modes m = 
0, 1, 2 for uj p /2tv = 10 GHz (c) and u p = (d). Cj = 0.15 pF, 
Co = 0.16pF/mm, D — 5nm, and W = 10 /im. 



which in this case reduces to 
2Cj 



{^1 -uj 2 m ) cot\k m e\ = k n 



(17) 



and the constant A rn (whose explicit expression is not 
given) is determined by the normalization condition (7). 
The coupling strength g m in Eq. (15) takes the form 



9n 



m COs(&. m £) 

DW 



2w m Cr 



(18) 



Figure 2 shows u m and g m as a function of the length 
t for a few lowest modes. It demonstrates that the cou- 
pling strength is in the strong coupling regime (g /2n ~ 
10 MHz) when the resonator is close to resonance to the 
diamond NV center (uj m ss e). Notice that a capacitive 
junction gives coupling as strong as a Josephson junction. 
A submicron size of the gap junction has been made on 
coplanar waveguide [24]. 

Double NV center. — Let us now consider the case 
with two NV centers. For simplicity, we assume symmet- 
rically located NV centers with subresonators of lengths 
£,£',£ in this order (xs = — xq = £ + £'/2, X4 = —x 2 — 
£' /2, and D = x 2 — X\ = X4 — x 3 <C £, £')■ Because of the 
symmetry, a normal-mode wavefunction ijj m has either 
even or odd parity. Here, we only consider the lowest 
mode, i.e., the A/2 mode. The eigenfunction takes the 
form of 

\cosk(x-£-£'/2) (£'/2<x<£) 



The mode frequency u> 
conservation (6) 



vk is determined by the current 



vCqlo = Cj (uj 2 p - lj 2 ) [cot(k£) - tan(M'/2)] (20) 
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FIG. 3: The coupling strength (a) for uj p /2it = 10 GHz and 
4otai = 18 mm and (b) for Ld p /2n = and £ to tai = 25 mm. 
Cj = 0.15 pF, Co = 0.16pF/mm, D = 5nm, and W = 10 /un. 



while the constants A (not given here) is determined from 
the normalization condition (7). The coupling strength 
at the junctions is found to be 
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WD V 2C E w 



[cos(fcf) - sm(k£) tan(fcf/2)] (21) 



The coupling strength g as a function of £' is shown in 
Fig. (3). 

Let us now estimate the coupling strength J between 
two NV centers through a virtual excitation of the pho- 
ton in a dispersive regime. For the resonators coupled 
through Josephson junctions, we have the frequency, 
uj/2tt — 2.4 GHz, and the NV center-resonator cou- 
pling, g/2n ~ 8 MHz, for £ total = 11 + 1' = 18 mm and 
£' ~ 0.4 mm. Assuming a detuning of A/27T = 200 MHz 
between e and u, we obtain the coupling strength be- 
tween the NV centers J = g 2 /A ~ 2tt x 320 kHz. For 
Q = 10 5 , the cavity decay rate is k = uj/Q ~ 2ir x 20 kHz 
which is lower than the coupling strength J. Even for 
J ~ k, two-qubit quantum gates can be performed with 
a high-fidelity as demonstrated in Ref. [7]. Thus, given 
the spin coherence time of 1 — 10 ms [25], a high fidelity 
quantum gate is realizable using our scheme. 

Discussion - Above we have focused on how to strongly 
couple diamond NV centers to the resonator. The strong 
coupling opens another interesting possibility to couple 
NV centers to other types of superconducting qubits [6, 
8-10, 16, 17, 26] through the resonator. This way one 
can take the best of features that each qubit provides. It 
will ultimately provide a novel architecture for quantum 
information processors integrating diamond NV centers 
and superconducting qubits into the circuit QED system. 
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